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We propose a method to decouple the nanomechanical resonator in optomechanical systems from 
the environmental noise by introducing a chaotic coherent feedback loop. We find that the chaotic 
controller in the feedback loop can modulate the dynamics of the controlled optomechanical sys¬ 
tem and induce a broadband response of the mechanical mode. This broadband response of the 
mechanical mode will cut off the coupling between the mechanical mode and the environment and 
thus suppress the environmental noise of the mechanical modes. As an application, we use the 
protected optomechanical system to act as a quantum memory. It’s shown that the noise-decoupled 
optomechanical quantum memory is efficient for storing information transferred from coherent or 
squeezed light. 

PACS numbers: 03.67.Pp, 02.30.Yy 


I. INTRODUCTION 

Optomechanical systems have attracted intense atten¬ 
tion in recent years due to its extensive applications 
and rapid progress has been made both theoretically and 
experimentally in related fields dH. One of the most 
interesting problems for optomechanical systems is to ex¬ 
plore the quantum aspects of mechanical motion [l^ - [l5l| , 
which is important not only for fundamental studies of 
quantum mechanics, but also for further applications, 
such as the detection of gravitational waves [ii[i3 , and 
quantum memorise [l8l.ll9l[. However, these quantum ef¬ 
fects will be damaged by environmental noises. Although 
the recent development of experimental techniques have 
made it possible to cool mechanical modes to the ground 
state [l2, [20l - l^ , the mechanical quantum superposition 
state is still too fragile under environmental noises, and 
thermal noise will be dominant if the mechanical mode 
is far from the ground state. 

Due the problems mentioned above, how to suppress 
the environmental noises more efficiently is crucial in ex¬ 
ploring the quantum-classical boundary of nanomechani¬ 
cal resonators. One possible way to solve this problem is 
to introduce either active or passive feedback to compen¬ 
sate the noise effects [20l - l^ . Side band cooling [20l - [25j 
is the most widely-used passive compensation method, 
and experiments |20l - l23l| in both the strong and the weak 
optomechanical coupling regimes have been rep orted to 
reach the quantum-mechanical ground state |20l - l23ll. Ap - 
proaches based on active feedback compensation (261 - 1^ , 
are also effective in suppressing environmental noise. The 


*Electronic address: jing-zhang@mail.tsinghua.edu.cn 


essence of these methods is to steer the system to the de¬ 
sired state by using the measurement output from a par¬ 
ticular quantum nondemolition measurement. Another 
possible way to solve this problem is to decouple the 
mechanical resonator from the heat bath by introducing 
a carefully-designed open-loop control. Dynamic decou- 
pfing control (DDC) and its optimized versions HH- 
l33 | are possible ways to achieve this, which introduce 
high frequency control pulses to average out the low fre¬ 
quency noises. However, it is not easy to generate the 
required high-frequency or optimized pulse in optome¬ 
chanical systems and, thus, to our knowledge, DDC has 
never been used to protect the mechanical states in such 
systems. 

Motivated by the DDC-type control and especially our 
recent work |35l| (introducing a broadband chaotic control 
to suppress decoherence of a superconducting qubit [H), 
in this paper, we propose a method to decouple the 
nanomechanical resonator from its environmental noises 
by introducing a chaotic coherent feedback loop. Based 
on the theory of coherent feedback [Sbl - ld^ . which is one 
of the major quantum feedback approaches [ttI - Is^ I. the 
basic idea of our method is to transfer a broadband 
chaotic control signal from the controller to the controlled 
optomechanical systems by feedback connections. This 
broadband control induces an effective broadband fre¬ 
quency shift of the mechanical resonator and then decou¬ 
ples the mechanical mode from the environmental noises. 
Afterwards, we use the protected mechanical mode as 
a quantum memory to store continuous-variable quan¬ 
tum signals, such as coherent states and squeezed states, 
which may have potential applications. 

This paper is organized as follows. In Sec. H, we will 
provide general discussions to show the noise-decoupling 
mechanism for our chaotic feedback strategy. The pos- 
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sible physical implementations for our noise-decoupling 
strategy in on-chip optomechanical systems are discussed 
in Sec. III. As an application, in Sec. IV, we show how to 
use a optomechanical system, protected by the designed 
chaotic feedback control, to act as a quantum memory. 
In Sec. V, we summarize the conclusions and provide a 
few forecasts of future work. 


II. NOISE DECOUPLING BY CHAOTIC 
FEEDBACK 


In this section, we show the mechanism of our chaotic- 
feedback-induced noise decoupling strategy, in particular 
for quadratically-coupled optomechanical systems HH- 
[^ . This is motivated by our previous work [35j | which 
shows that decoherence in supercoducting circuits can be 
greatly suppressed by chaos which is typically believed 
to be a source of decoherence. The main idea of the 
chaos-induced decoherence suppression approach is to in¬ 
troduce a broadband chaotic signal to ’’randomly” kick 
the system and compensate the effects of noise. This 
is somewhat similar to the noise suppression approaches 
by the quantum Zeno effect in which random signals are 
introduced to kick the system to compensate the noise 
effect. However, chaotic signals are deterministic signals 
and thus will not introduce additional decoherence. 

Note that, there are some difficulties in introducing 
such kind of chaotic control to suppress the noises of the 
quantum-mechanical mode in optomechanical systems: 
(i) it is quite hard to drive the mechanical mode of an op¬ 
tomechanical system directly by a chaotic acoustic field; 
and (ii) the optical cavity in the optomechanical system 
will work as a low-pass filter to squeeze the broadband 
chaotic signal if we drive the system directly by an open- 
loop chaotic optical signal and thus make the control sig¬ 
nal not so ” random”, which would lead to a failure of our 
decoherence-suppression approach. To solve these prob¬ 
lems, we introduce a particular coherent feedback loop to 
break the symmetry of the optomechanical system. Thus, 
the chaotic controller in the feedback loop can broaden 
the bandwidth and preserve the high-frequency compo¬ 
nents of the mechanical mode, and protect it from the 
environmental noises. 


As illustrated in Fig. [H our feedback control system 
consists of two components, i.e., a quadratically-coupled 
optomechanical system (the controlled system) and a 
chaotic controller. These two components are connected 
by a mediating optical field, from which we can con¬ 
struct a field-mediated coherent feedback system [3^ 
IHi 13 ■ In the interaction picture of the noise fre¬ 
quency w, the Hamiltonian of the controlled device, i.e., 
the quadratically-coupled optomechanical device, can be 
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FIG. 1: (Color online) Schematic diagram of the noise¬ 

decoupling system by coherent feedback modulation. Two 
quantum components, i.e., a quadratically-coupled optome¬ 
chanical system and a chaotic controller, are connected by 
the mediated optical fields. The output of the optomechani¬ 
cal system is taken as the input fed into the chaotic controller. 
Also, the chaotic signal generated by the chaotic controller is 
then fed back to control the dynamics of the quadratically- 
coupled optomechanical system. The Ml, M2, M3, M4 repre¬ 
sent total-reflection mirrors that are introduced to change the 
light path. In this model, the radiation pressure can directly 
change the frequency of the mechanical resonator because the 
coupling between the mechanical mode and the optical mode 
is quadratic. 


written as [Sll - I^ 
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where oi and bi denote the annihilation operators 
of the cavity mode and the mechanical mode in the 
quadratically-coupled optomechanical system, and 
Hi are the natural frequencies of these two modes. Here, 
we assume that h = 1. The optomechanical coupling we 
consider here is a kind of quadratic optomechanical inter¬ 
action with strength Gi The optical mode oi is 

driven by an external driving field with strength £i and 
frequency . Here b{uj) represents the noise mode with 
frequency uj acting on the mechanical mode and g{uj) is 
the coupling strength between the mechanical mode and 
the noise mode. 

Here we use to denote the Hamiltonian of the 
chaotic controller, and 02 denotes the annihilation op¬ 
erator of the chaotic cavity field in the controller. Then 
the interaction Hamiltonian of the quadratically-coupled 
system and the controller Hint takes the form (see Ap¬ 
pendix A) 


Hint = ^(\/7i 72 - \/727/)(4ai “ 4«2), 


( 2 ) 
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where 71 and 72 represent the damping rates of the opti¬ 
cal cavities in the controlled system ” 1 ” and the chaotic 
controller ” 2 ”, and 7 / denotes the damping rate of the 
controlled cavity induced by the feedback field. The total 
Hamiltonian of the coherent feedback loop is provided by 


which the mechanical mode is efficiently decoupled from 
the environmental noises. 

As shown in Appendix B, the correction factor M is 
determined by the power spectrum Sf^uj) of the chaotic 
signal /(t) 


Htot = Hi+H, + Hint. (3) 


In the strong-driving regime, the optical fields in the 
quadratically-coupled optomechanical system and the 
chaotic controller can be treated classically. Here we re¬ 
place the operator oi by ai(t), which represents the clas¬ 
sical part of the optical field oi, and then eliminate the 
classical parts including He and Hint in tbe total Hamil¬ 
tonian. Thus the Hamiltonian of the feedback control 
system given in Eq. © can be simplified as 


Heff =flibib\ + f{t) b\hi 

+ '^g{uj) b{uj)b\ exp{iujt) +'h.c. 


(4) 


where f{t) = G'i|ai(t)|^, and the amplitude of the cavity 
field |ai(t)| is modulated by the chaotic controller and 
thus it is a broad-band signal. The effective Hamiltonian 
in Eq. (|4|) includes three parts: (i) the free Hamiltonian 
of the mechanical mode with natural frequency Hi; (ii) a 
correction term with the mechanical frequency shift f{t) 
induced by the chaotic controller He] (hi) the interaction 
Hamiltonian Hint between the mechanical mode 61 and 
its environmental noises b{ui). In the rotating reference 
frame with the unitary operator 


U = exp 



(/(r) -b Hi) b\bidT 


(5) 


the effective Hamiltonian is given by 


Hes = H^HesU -iU'fdU/dt 


g{oj) + h.c. 


( 6 ) 


By averaging over the broadband signal f{t) we have 
(see Appendix B) 


exp 


-i / f{T)dT 
^0 


= \/m, 


(7) 


where M is a correction factor. Thus, the effective Hamil¬ 
tonian shown in Eq. © can be simplified as 


HeS = 


g(w) | 6 (w) 6 | exp[—i(Hi — uj)t] + h.c.j , ( 8 ) 


where g{ui) = \fMg{u}) is the modified coupling strength 
between the mechanical mode and the heat bath after 
introducing the chaotic signal f{t). It can be seen that 
the modified coupling strength g(u}) can be greatly de¬ 
creased if the correction factor M is small enough, under 


M = exp 



/,,2 



(9) 


where uJu and w/ are the upper bound and lower bound of 
the frequency band of the chaotic signal / (t). Note that, 
M varies from 0 to 1. Specially, M = 0 corresponds 
to the full-decoupling case, and M = 1 corresponds to 
the case without decoupling. Since the power spectrum 
Sf{uj) is broadened by the chaotic modulation, the value 
of M is thus very small and the mechanical mode is de¬ 
coupled from the environmental noises. 


III. PHYSICAL IMPLEMENTATION IN 
ON-CHIP OPTOMECHANICAL SYSTEMS 


In this section, we discuss how to physically implement 
our chaotic-feedback-based noise decoupling strategy in 
on-chip optomechanical systems. 


A. Implementation of the quadratically-coupled 
optomechanical system 


Here we list two possible exaniples of the quadratic¬ 
coupling optomechanical system [5ll - l55j . The first exam¬ 
ple is shown in Fig. [2Ja), in which a membrane is placed 
in the middle of a cavity and can move freely under the 
laser-induced pressure [5ll454l |. Such kind of structure 
leads to a quadratic coupling term between the mechani¬ 
cal mode and the cavity mode. Another example for the 
quadratic-coupling is the rectangular membrane optome¬ 
chanical system [^. As seen in Figs, ^b) and (c), the 
rectangular membrane placed above a toroidal cavity is 
driven by the optical field inside the toroidal cavity, which 
may generate both linear coupling and quadratic cou¬ 
pling modes between the cavity field and the membrane. 
The coupling strengthes of these two coupling modes are 
determined by three factors: (i) the vibrational mode 
of the rectangular membrane; (ii) the distance between 
the membrane and the upper surface of the toroidal cav¬ 
ity; and (iii) the relative position of the toroidal cavity. 
Moreover, the coupling modes displayed in the rectangu¬ 
lar membrane optomechanical system can be controlled 
by modulating the above factors. The purely quadratic¬ 
coupling mode can be realized when [^: (i) the rect¬ 
angular membrane is excited in a vibrational mode that 
contains at least one node; (ii) the rectangular membrane 
is placed right above the toroidal cavity; and (iii) the 
node of the membrane is located at the central point of 
the cavity. Under these conditions, the linear coupling 
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FIG. 2: (Color online) Schematic diagram of the quadratic¬ 
coupling optomechanical systems with a Fabry-Perot cavity 
and the rectangular membrane, (a) Quadratic optomechani¬ 
cal system with a Fabry-Perot cavity: the quadratic-coupling 
is realized by putting a membrane in the middle of the Fabry- 
Perot cavity, (b) Top view and (c) cross-sectional view of a 
rectangular membrane optomechanical system, where its node 
coincides with the central point of the cavity. The rectangu¬ 
lar membrane supports various vibrational modes u = {j,k), 
where j,k = 1,2... are the mode indexes. Here the rectan¬ 
gular membrane is driven to the ( 1 , 2 ) mode, which has two 
anti-nodes and one node. 


term between the membrane and the cavity field can be 
completely removed. 

The mechanism of the rectangular membrane op¬ 
tomechanical system is similar to the Fabry-Perot-type 
quadratic-coupling system, and they share the same 
Hamiltonian, which is shown in Eq. ©• Hereafter, we 
apply our noise-decoupling method to the rectangular 
membrane optomechanical system presented above. 
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FIG. 3: (Color online) The noise-decoupling model with 
the control of a toroidal cavity. Here the toroidal cavity 
is a chaotic controller, which shifts the cavity field of the 
quadratically-coupled optomechanical system to chaos. 


B. Implementation of the chaotic controller 


In this section, we consider an optomechanical system 
[see Fig. [3] with chaotic dynamics (s^j as the chaotic 
controller in the feedback control loop. For simplicity 
we denote the controlled quadratically-coupled optome¬ 
chanical device as system 1 , and the chaotic controller as 
system 2. The Hamiltonian of system 1 is displayed in 
Eq. (IT}; and the Hamiltonian of system 2 is taken as 

H2 =^02^2^2 + G2 0202(^2 “f ^2) + f^2 ^2^2 (10) 

-b ie2[a\ exp{—iuJd2t) — 02 exp{iujd2t)], 


where a 2 and 62 denote the annihilation operators of the 
cavity mode and the mechanical mode in system 2 ; and 
Was I correspond to their inherent frequencies. Here G 2 
denotes the optomechanical coupling strength in system 
2. The cavity mode in system 2 is driven by an input 
laser field with driving strength £2 and corresponding 
driving frequency uJd 2 - Here, the driving frequencies of 
the cavity modes in the two systems are chosen to be: 
ujdi = 1-0d 2 = 00 d- In the rotating reference frame with 
the unitary operator U = exp[—iujd{a\ai + a\a 2 )t], the 
total Hamiltonian of the quantum feedback loop can be 
transformed to the form 

Htot =Aia|ai -|- Gia\aib\bi Q16J&1 

-|- A2Cl2a2 “b G'2a2®2(^2 “f ^ 2 ) “b 142^2^2 
+ i£i{a\ — ai) + i£ 2 {al — a 2 ) 

1 4 . 4 . ■^/ 

+ ^(\/7i72 - \/727/)(4 «i - 4«2) 

+ g{Lo)[b'^{Lo)bi exp(—fwt) -b h.c.]. 
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where Ai = uja^ — ujd, and A 2 = — cud, denote the 

detuning frequencies of cavities 1 and 2. Here 71 and 
72 represent the damping rates of the optical cavities 
1 and 2 , 7/ denotes the damping rate induced by the 
feedback field of the controlled cavity. We use the quan¬ 
tum Langevin equations to describe the dynamics of the 
chaotic feedback system 

hi = - lAiOi - “ iGiaib\bi 

- ^lllf 0.2 + £1 - (\/q 7 + V^)ai 4 n, 

02 = — —^02 — iG 2 a 2 (bl -|- 61) -I- £2 

— \/7l72 Ol — 


( 12 a) 

( 12 b) 


Yi 

bi = -i^libi - iGia\aibi - —bi - \/l\ (12c) 

62 = —i^ 2^2 ~ iG 2 C^ 2 ^ 2 , -^^2 ~ '\/r^^2,in) (12d) 

where oi^n ( 02 ,in) is the input of the optical cavity in 
system 1 ( 2 ); 61 ,in (^ 2 ,in) and Fi (r 2 ) are the input and 
the damping rate of the mechanical mode in system 1 ( 2 ). 
We assume that the backaction of the mechanical mode 
acting on the optical mode in system 1 is very weak, 
then the evolution of the cavity mode 1 mainly depends 
on Eqs. (Hi), dUb), and (IT^ l. In the strong-driving 
regime, the semiclassical approximation can be applied: 
oi = ai -|- oi, 02 = 0:2 + 02 , and 62 = /52 + & 2 , where aj, 
02 , and /32 represent the classical parts and oi, 02 and 62 
denote the operators for the quantum fluctuations. Then 
we neglect the quantum fluctuation terms in Eqs. (I12h '). 
(H^), and (I12h l . Thus the evolution of the classical parts 
in the total system can be described by 

di =-iAiai- 

£1 — ^727/ 0:2, 


d2 = — 1^202 — ^Ci2 — ^G20L2(fi2 + ^2) 
+ £2 - VTi^ai, 


(13b) 


/32 


—iniP2 — iG2a2a2 



(13c) 
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FIG. 4: (Color online) Power spectra of the mechanical mode 
in system 1: (a) without feedback and (b) with the chaotic 
feedback is applied, (a) the mechanical membrane is driven 
by a period signal; (b) we use a chaotic controller (optome¬ 
chanical system in this case) to modulate the power of the 
mechanical membrane to a wide region. The parameters are 
set as follows: Ai/27r = 0.75 GHz, A2/27r = 0.12 GHz, 
7i/27r = 1 MHz, 72/27r = 0.24 GHz, 7//27r = 0.05 MHz, 
ri/27r = 0.01 MHz, T 2 = /27r = 1.4 MHz, !2i/27r = 1 MHz, 
n2/2TT = 0.345 GHz, G 1 / 2 -JT = 0.1 MHz, G2/2tt = 0.1 MHz, 
£i/27r = 6.6 GHz, and £2/27r = 13.2 GHz. 


When the strength of the driving field £2 is strong 
enough, the optomechanical system 2 enters the chaotic 
regime and will have a broadband cavity spectrum. As 
the chaotic controller, system 2 can spread the spectrum 
of system 1 both in the cavity mode and in the mechanical 
mode. Figure. U shows the spectrum of the mechanical 
mode in system 1 without [Fig.UKa)] and with [Fig.|TKb)] 


the feedback modulation. As shown in Fig. Idja), only a 
single peak with very small sidebands is displayed in the 
spectrum of the mechanical mode if we do not introduce 
any feedback modulation. The power of the background 
frequency components is very small (less than -150 dBm). 
This corresponds to the periodic case. After we intro¬ 
duce the chaotic feedback [see Fig. SKb)], the spectrum 
of the controlled mechanical mode is greatly broadened 
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and the whole baseline of the spectrum is increased to 
above 150 dBm. This corresponds to the chaotic case, 
and the broadband response of the mechanical mode will 
decouple the mechanical mode from the environmental 
noises. 

As discussed in Sec. II, we use the factor M to eval¬ 
uate the efficiency of our noise decoupling strategy [see 
Eq. (ini)]. The value of M is determined by the spectrum 
Sf{io) of the signal f{t) (recall that f{t) = Gi|ai(t)|^), 
which can be obtained by numerically solving Eq. (USD. 
Note that M ^ 1 when the spectrum Sf{uj) is concen¬ 
trated in a narrow region, and M will be close to zero if 
the spectrum Sf (aj) is broadened by the chaotic modula¬ 
tion. In our numerical simulations, we find that M k, lil 
we do not introduce feedback [Fig.SKa)] and M = 0.0074 
if we introduce the chaotic feedback [Fig. |ll[b)], which 
coincides with what we expect. 

IV. STORAGE OF CONTINUOUS-VARIABLE 
QUANTUM INFORMATION 

The storage of continuous-variable quantum informa¬ 
tion, i.e., to realize continuous-variable quantum mem¬ 
ory [ll, [I^, [5§ - l^ . is important for quantum com¬ 
munications and quantum computation. One possible 
way to solve this problem is to transfer the continuous- 
variable information in the optical signal to an on-chip 
mechanical resonator which has a lower damping rate. 
The continuous-variable optomechanical quantum mem¬ 
ory system we consider here is presented in Fig. |5l which 
includes the input (output) fields, an optical cavity, and 
a mechanical resonator [1^. By exchanging states, be¬ 
tween the cavity mode and the mechanical mode, a quan¬ 
tum state carried by the input held can be written into 
and stored in the nanomechanical resonator. 

However, the quantum information stored in the me¬ 
chanical resonator will unavoidably be destroyed due 
to the coupling between the mechanical resonator and 
the environmental noise. Thus, to realize such kind of 
continuous-variable quantum memory, we have to sup¬ 
press the decoherence effects of the mechanical mode 
induced by the environmental noise. As we have dis¬ 
cussed in the previous sections, introducing a chaotic co¬ 
herent feedback loop to drive the mechanical mode into 
the broad-band regime is an efficient way to decouple 
the mechanical mode from the environmental noise. In 
this section, we will show how to use this noise-decoupled 
nano-mechanical resonator as a quantum memory. 


Our purpose here is to use a noise-decoupled mechan¬ 
ical resonator to store continuous-variable information. 
The key point is how to transfer a quantum state to a 
mechanical mode and decouple this mechanical mode si¬ 
multaneously. Here we propose a strategy with two op¬ 
tical cavities sharing the same mechanical resonator but 
with different optomechanical coupling: one is with linear 
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FIG. 5: (Color online) Schematic diagram of an optmechani- 
cal system for quantum information transfer and storage. A 
beam of light with a desirable quantum state is fed into a cav¬ 
ity, and then transferred to the mechanical resonator. Here 
a is the cavity mode and b denotes the mechanical mode, n 
represents the mean thermal excitation phonon number which 
follows the Boltzmann distribution. 

optomechanical coupling used for quantum memory; and 
the other is with a quadratic optomechanical coupling, 
used for noise decoupling. 


Let us now consider how to apply this quantum mem¬ 
ory model in the rectangular membrane optomechani¬ 
cal system proposed in Ref. [s^. Figure ini[a) shows two 
toroidal cavities (A and B) connected to a rectangular 
membrane. The types of coupling between the cavity 
mode and the mechanical mode are determined by the 
positions they are placed: the node of the membrane 
corresponds to a linear coupling and the anti-node cor¬ 
responds to a quadratic coupling [Fig. UKa)]. Thus, we 
place the toroidal cavity (cavity A) used for noise de¬ 
coupling at the node of the membrane; and the other 
toroidal cavity (cavity B), used for quantum memory, at 
the anti-node. Toroidal cavity A is modulated by the 
chaotic controller (toroidal cavity C), which leads to the 
decoupling between the membrane and its environmen¬ 
tal noises. The cavity B is used for storing the quantum 
state in the membrane. The coupling between the cavity 
mode and the mechanical mode is assumed to be lin¬ 
ear under the strong-driving regime [Tl,[l^. Thus, the 
Hamiltonian of the total system can be written as 

H =Asalas + Gsiasb\ + aj6i) 

-f (Hi -b Gi|ai(t)|^)6i^6i 
+ ^ g{uj)[b'^{u!)bi exp(—ftut) + h.c.j, 

iJ 

where Ug (aJ) represents the annihilation (creation) op¬ 
erator of the optical mode in cavity B, and LOg is the 
corresponding inherent frequency. Here Ag = ujg — ujd 
is the detuning frequency of cavity B, and uJd is the fre¬ 
quency of the external driving field. Also, Gg denotes 
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FIG. 6: (Color online) (a) Top view of the quantum mem¬ 
ory system. The noise-decoupled quantum memory system 
can be divided into two parts shown by the dashed rectan¬ 
gular grid frames: the chaotic feedback loop (inside the blue 
frame) for the noise decoupling of the rectangular membrane; 
the setup used for transferring the quantum state (the red 
frame) from the input light to the noised-decoupled rectan¬ 
gular membrane, (b) Cross-sectional side view of the rectan¬ 
gular membrane optomechanical system. Cavity A is placed 
at the node of the rectangular membrane, and cavity B is 
placed at the anti-node. 


the coupling strength between the optical mode and the 
mechanical mode. To compensate the effect induced by 
the chaotic feedback on the quantum memory system, we 
take the detuning frequency as Ag = fli + Gs\ai{t)\'^. In 
the rotating reference frame with the unitary matrix 


U = exp 



{Gi\ai{T)\'^ + ni){b\bi + a\ai)dT , 

(15) 


the effective system Hamiltonian can be represented by 
i7eff = Gg(al6i+ag6l)+^5(u;)[6+(u;)&ie-*(^^-‘")‘+h.c.], 

CO 

(16) 

where g{oj) = \/Mg{uj) and M is the decoupling factor. 
After introducing the adiabatic approximation to elim¬ 
inate the cavity mode shown in Ref. [l^ . we use bi to 
denote the annihilation operator of the mechanical mode 
and the quantum Langevin equation of the optomechan¬ 
ical system can be simplified as 

^ - y/vQd-\/l\bin{t), (17) 

Qd denotes the optical field fed into cavity B. Let ad = 
ctd+bdj where ad and ad denote the classical part and the 
quantum fluctuation of the optical mode. The fluctuation 
terms dd and 6in satisfy the relations: {dd{t)d\{t )) = 
{bin{t)bl^{t )) = {n + l)S{t-t), where n (Hi) Ri 
ksT/Mil is the mean thermal excitation phonon number. 
The parameter v in Eq. (II3 can be calculated by = 
(GslodD^/Ts, where Gg is the coupling strength between 
the mechanical mode and the optical mode, and 7 s is the 
damping rate of the optical mode [T^ . 

We now assume that the system is initially in a Gaus¬ 
sian state. We use the fidelity Foo between the initial 
state and the steady state of the mechanical mode to 
characterize the efficiency of noise decoupling, which can 
be calculated by [3 


=(«'o|Poo|«'o) 


n 

i=±s 


exp(j) 


ri(2n + 1- exp(j)) 


2{iy + Ti) 


(18) 


Here s is the squeezing factor (see Appendix C). The 
steady-state fidelity Foo mainly depends on four factors: 
the mean thermal excitation phonon number n, the cou¬ 
pling strength the squeezing factor s, and the mechan¬ 
ical damping rate R. We can see that the fidelity Foo can 
be increased by decreasing the mechanical damping rate 
Ti, and, as shown in Sec. H, R can be reduced by intro¬ 
ducing a chaotic feedback loop. In fact, after introducing 
the chaotic feedback control, the effective damping rate 
of the mechanical mode is given by 


r; = MTi. 


(19) 


Thus the modified fidelity F/^ can be written as 


F' 


n 


exp(j) -b 


r'^(2n + 1- exp(j)) 

2iu + r[) 


1 

2 


( 20 ) 


When the controller in the feedback loop enters the 
chaotic regime, we have T'^ r: 0 , and thus F/^ r 1 , which 
means almost perfect quantum state transfer. 

Then, we numerically calculate the steady-state fidelity 
Foo between the input state and the steady state of 
the mechanical resonator. Two different Gaussian in¬ 
put states are considered: coherent states and squeezed 
states. 




































A. Coherent input state 


In this subsection, we consider the quantum mechan¬ 
ical memory system with a coherent input state. For a 
coherent input state, the squeezing factor s = 0. Thus, 
in this case, the fidelity can be simplified as 





Fin 
V + 


-1 


( 21 ) 


By comparing the fidelity between the input state 
and the steady state of the mechanical mode (under the 
noise-decoupling control [see Fig. [7l[a)] and without the 
noise-decoupling control [see Fig. [3(b)]), we find remark¬ 
able improvement of the efficiency of the quantum mem¬ 
ory by introducing chaotic control. From Fig. |3[a) and 
Fig. 13(b), we can observe that the decrease of the mean 
thermal excitation phonon number n or the increase of 
the parameter u would lead to the improvement of the 
fidelity of the quantum transfer. If we fix the parameter 
zz = 50 kHz, the fidelity of the quantum transfer will fall 
to zero rapidly when increasing the excitation phonon 
number n without introducing the noise-decoupling con¬ 
trol [Fig. [3(a)]. We find that the fidelity of the quan¬ 
tum memory is increased and approaches one even when 
the mean thermal excitation phonon number n exceeds 
10^ after introducing the noise-decoupling control. This 
means that our noise-decoupling method efficiently re¬ 
duces the damping rate of the mechanical mode Fi, and 
thus protects the coherent input state from decoherence. 


B. Squeezed input state 


CD 

T3 



1 
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n/1000 
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CD 
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Let us consider the case that the input state is a 
squeezed state with squeezing factor s ^ 0. By adjusting 
the squeezing factor s and the mean thermal excitation 
phonon number n, we study the fidelity between the in¬ 
put squeezed state and the steady state of the mechanical 
mode. 

Compared to case without noise-decoupling control 
shown in Fig. [8](a), the fidelity under noise-decoupling 
control is significantly improved [see Fig. [HKb)] for differ¬ 
ent chosen system parameters. As shown in Fig.[8Ka) and 
(b), the fidelity decreases when increasing the squeezing 
factor s and the mean thermal excitation phonon num¬ 
ber n. Here we vary the squeezing factor s from —5 to 
5, and it can be found that the curve of fidelity is sym¬ 
metrical about the plane s = 0 in the three-dimensional 
fidelity space. For each parameter n, the fidelity is max¬ 
imized when s = 0, which corresponds to the case that 
the input state is a coherent state. The quantum infor¬ 
mation stored in the memory system is more likely to be 
damaged by the heat bath when increasing the degree of 
the squeezing factor s. As shown in Fig. [3 the fidelity of 


FIG. 7; (Color online) The fidelity (a) before the noise de¬ 
coupling and (b) after the noise decoupling. Here n is the 
mean thermal excitation phonon number which follows the 
Boltzmann distribution, and is a parameter related to 
the optomechanical coupling strength. The parameters are: 
Hi/27r = 1 MHz, ri/27r = 5 Hz for (a), and r'i/27r = 0.037 
Hz for (b). 


quantum transfer F = 0.16 is very low when n = 10® and 
s = 0 without the noise-decoupling control [see Fig.[8Ka)], 
while, with the same condition, the fidelity is enhanced to 
he F = 0.96 if we introduce the noise-decoupling control 
[see Fig. [3(b)]. When the squeezing factor s is increased 
to approach 5, the fidelity decreases to zero rapidly with¬ 
out the noise-decoupling control [see Fig. [3(a)], while it 
will remain nonzero, i.e., F = 0.38, when we introduce 
the noise-decoupling control [see Fig. [3b)]. 
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n/1000 Squeezing factor s 



0 . 2 - 


n/1000 


Squeezing factor s 


5 


FIG. 8: (Color online) The fidelity (a) before noise decoupling 
and (b) after noise decoupling. The parameters used are: 
i//27r = 10 kHz, ri/27r = 5 Hz for (a), and r'i/27r = 0.037 Hz 
for (b). The natural frequency of the mechanical mode is 
assumed as fli/27r = 1 MHz. 


V. CONCLUSION 


To summarize, by introducing a chaotic feedback con¬ 
trol loop, we propose a strategy to decouple a nanome¬ 
chanical resonator in a quadratically-coupled optome¬ 
chanical system from the environmental noises. The 
main advantage of this method is to introduce a chaotic 
controller to significantly broaden the spectrum of a me¬ 
chanical resonator and thus efficiently suppress the envi¬ 
ronmental noise. As an application, we study this pro¬ 
posed the noise-decoupled nanomechanical resonator un¬ 
der chaotic coherent feedback control as a quantum mem¬ 
ory to store the information transferred from external 
optical signals. Two different input states, i.e., coherent 
and squeezed states, are studied to show the efficiency of 


the quantum memory. The numerical results show that 
the fidelity of the quantum memory has been greatly im¬ 
proved after introducing our noise-decoupling strategy. 
We believe that this nonlinear coherent feedback strategy 
will have various applications, such as nonlinear modu¬ 
lation of photon transport and high-sensitivity quantum 
measurements, which will be considered in future work. 
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Appendix A: Theory of Markovian coherent 
feedback network 

To study the multi-channel quantum input-output net- 
workjwe now introduce the SLH method presented in 
Ref. [^. In the SLH language, an open quantum sys¬ 
tem can be fully characterized by G = (S', L, H), where S 
denotes a nx n unitary scattering matrix, which satisfies 
SS^ = = I, L represents the dissipation operator 

which is determined by the dissipation channels induced 
by the input fields, and H is the free Hamiltonian of the 
system. Within the framework of G = {S,L,H), the 
quantum Langevin equation of an arbitrary system op¬ 
erator X is given by 

A = - + {Lt[A,L] + [L\X]L}/2 

+ {k4L\X] + [X,L]bl}. 

The SLH method provides a convenient way to study 
the all-optical quantum coherent forward and feedback 
networks [^ . Eor example, we show in Pig. two 
quantum components: Gi = {Si, Li, Hi) and G 2 = 
{S 2 ,L 2 , H^). The series product of these two components 
can be parameterized by 
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G2 > Gi — [<S'2iS'i, L2 + <S'2ii, 

Hi + H 2 + —{l\S2Li — l|S'|L 2)]- (A2) 


Sub-system 1 Sub-system 2 



FIG. 9: (Color online) Schematic diagram of the series prod¬ 
uct of two cascaded-connected components. 


Sub-system 1 



A typical coherent feedback control system is shown in 
Fig.Hni which is composed of the controlled system, i.e., 
system 1, and the controller, i.e., system 2. This coherent 
feedback control system can be seen as a series product of 
three components: Gi = (5i,Li,iJi), G2 = (52, T2, ^^2), 
and Gf = {Sf,Lf,Hi). Thus, the corresponding SLH 
parameters of this feedback control system can be repre¬ 
sented by 

G/> G 2 > Gi = {S,L,Hsys), (A3) 

where 

S = SfS 2 Si, L = S 2 S 1 L 1 S 1 L 2 + Lf, (A4a) 


FIG. 10: (Color online) Schematic diagram of a coherent 
feedback loop. 

Accordingly, the quantum Langevin equations of the 
two cavity modes oi and 02 can be represented by 

ai=-i[ai,Hi +H2]-]^{y/^+y/^fai 

- y/l2lf 02 - (Vbh + Vh) «l,m, 

12 

0,2 = -i[a2, Hi + H2] - —a2 - ^7172 ai - \/^ 02 ,in, 

(A8b) 


H,y, = Hi+H2 + Hi^u (A4b) 

and the interaction Hamiltonian induced by the coherent 
feedback loop is given by 

Hint =^XLIS2 Li - l\sIL2 
1% 

-f L)SfL2 - LlsjLf + Lt5/52Li - L\slsjLf). 

(A5) 

As an example, let us consider our feedback-induced 
noise-decoupling system. As introduced in section 3, a 
quadratically-coupled optomechanical device (system 1) 
and a chaotic controller (system 2) are connected by op¬ 
tical fields to construct a coherent feedback loop, which 
is similar to that given by Eq. (HU. Let oi (02) be the 
annihilation operator of the cavity mode in quantum sys¬ 
tem 1 (2) with corresponding damping rate 71 (72), and 
7/ is the damping rate of the controlled cavity induced 
by the feedback field. In this case, we have Li = ai, 
L 2 = V^ ffl2, and Lf = and Si = S 2 = Sf = I. 

From Eq. (ESI), the dissipation operator of the total feed¬ 
back loop can be written as 

A = (\/^+v^)“i + (A6) 

and the total Hamiltonian of the quantum feedback loop 
can be obtained from Eq. (IA4I) and Eq. (IA5I) as 

Hsys = Hi + H 2 + iLint 

= Hi+H2 + ^(V7i72 - V727/)(4«i - aia2)- 

(A7) 


where oi^n (02,in) is the input field fed into the system 
1 (2). By substituting Eqs. (P) and (ITUl) into Eq. (IA8I) . 
we can obtain the quantum Langevin equation given by 

Eq. dni). 


Appendix B: Derivation of the decoupling coefficient 

M 


The decoupling coefficient M is determined by the clas¬ 
sical cavity field f{t), which can be decomposed into 
a series of frequency components by the Fourier trans¬ 
form [sslls^. 

00 

f{t) = '^AnCOs{uJnt + lfn), (Bl) 

n—0 

where ujn, A„ and denote the frequency, the amplitude 
and the initial phase of the n-th frequency components. 
Integrating /(t) gives the control-induced phase shift 


ft “ /i 

= / /('^) d-T = ^— sin(a;„t -f ipn). (B2) 
•^0 ^0 

By introducing the Bessel-series expansion, we have 


exp [—i9{t)] = exp 


A, 


-f — sin(w„t -I- ( fin ) 


n=0 


= ^ ^ j exp[-maa;at - iUaiPa], 

(B3) 
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where Jna is the n-th Bessel function of the first kind. We 
then neglect the high-order terms in Bessel series, which 
can be considered as the fast variables in the system, and 
only keeps the zero-order terms in Eq. (IB3I) , by which we 
have 


parameter. Here M and N satisfy the inequality > 
N{N -I- 1). Then we introduce the squeezing factor s 
of the input quantum state, which is given by 

s = \n[M + M* + 2N + 1]. (C5) 


exp{-i9{t)) = 


= exp 


y^ln Jq 



Under the condition that Aa Wq, the zero-order Bessel 
term can be approximately expressed as Jo{Aa/u!a) « 
1 — {Aa/2uJa)'^. Furthermore, from A^ oJa, we have 
ln(l — {Aa/2uJaY) Ri —{Aal2ujaY. Thus Eq. (IB4I) can 
be simplified as 


n^o 



exp 


exp 


1 



,.,2 



(B5) 


Let \/M = exp(—id(t)), and M is defined as the decou¬ 
pling factor, then from Eq. (IB5I) we have 


To calculate the fidelity of the quantum memory, let 
us define the normalized position a; = (6i -I- h\)/\/2, mo¬ 
mentum p = (bi — b\)/'/2i, and the conjugate vector 
z = {x,p) of the mechanical mode. We also introduce 
the symmetrized covariance matrix V, which is given by 

V = (C6) 

where Az = z — (z). With Ito’s rule d(ab) = (da)b + 
a(db) + da db, the time evolution of the covariance matrix 
V is described by the Lyapunov differential equation 

V = AV + VA'^ + ri{n+l/2)l2 + vA, (C7) 

where A = — [(^-|-ri)/2]/2, and I 2 is the two- 

dimensional identity matrix. Here, A is a matrix related 
to the degree of squeezing, which can be calculated by 


M = exp 



(B6) 


Appendix C: fidelity of the quantum memory 

The Langevin equation of the mechanical operator bi 
is shown in Eq. (113. The steady value of the mechanical 
mode can be obtained by setting dhi/dt = 0 as 

(61(00)) = —(Cl) 

1 / + ii 


I f 2N+1 + M + M* M - M* \ 

^~2\ M-M* 2N+1-{M + M*) ) ■ 

. 

For a squeezed input state, the fidelity between the initial 
state and the steady state of the mechanical mode is given 
by 


Foo =(«'o|Pcx.|«'o) = 


n 

j = ±S 


exp (j) 


y/ det(Uoo + Vo) 

ri(2n -I- 1 - exp(j)) 

2 (^ + ri) 


(C9) 


where (•) is the average over the input vacuum fluctu¬ 
ation. We then define the quantum Wiener processes 
A{t) = fg dd{t)dt , B{t) = fg bin(t')dt', by which we can 
obtain the quantum stochastic differential equation from 
Eq. (flTll as 

dbi = - —-bidt — y/v addt —y/v dA—y/VidB. (C2) 

The quantum fluctuation terms dA and dB satisfy that 

(cM) = {dB) = 0, (C3) 

and obey the quantum Ito rules 

dA dA"! = {N + l)dt, dA^dA = N dt, 

(dA)^ =Mdt, ((Mi)2 = MUt, (C4) 

dB dB^ = {n + \)dt, dB^dB = n dt, 

where n represents the thermal excition number, N is the 
effective photon number, and M denotes the squeezing 


where Uoo denotes the stationary solution of the Lya¬ 
punov differential equation [Eq. (IC7|) ] and Vb is the co- 
variance matrix of the input state, which can be calcu¬ 
lated by 


When the input state is a coherent state, such that M = 
N = 0 and thus s = 0, the fidelity in this case can be 
simplified as 




1 

y/ det(Uoo + Uo) 



Tin 
V + Ti 


-1 


(Cll) 


It can be found from Eqs. (jC9|) and (|C11|1 that the fi¬ 
delity increases when increasing the mechanical damping 
rate Ti for both squeezed states and coherent states. 
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